Magnetic-field dependence of electron spin relaxation in n-type semiconductors 
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We present a theoretical investigation of the magnetic field dependence of the longitudinal (Ti) and 
transverse (T2) spin relaxation times of conduction band electrons in n-type III-V semiconductors. 
In particular, we find that the interplay between the Dyakonov-Perel process and an additional 
spin relaxation channel, which originates from the electron wave vector dependence of the electron 
g-factor, yields a maximal T2 at a finite magnetic field. We compare our results with existing 
experimental data on n-type GaAs and make specific additional predictions for the magnetic field 
dependence of electron spin lifetimes. 
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The electron spin in a semiconductor is a robust ob- 
ject which can be utilized to add new functionality to 
existing electronic devices or to even build completely 
new devices based on this spin degree of freedom. Q Es- 
tablishment of successful spintronics devices requires a 
thorough understanding of the electron spin dynamics in 
a semiconducting environment. In particular, spin relax- 
ation processes need to be identified and controlled. 

Important electron spin relaxation processes in n- 
type semiconductors include the Elliott- Yafet (EY) pro- 
cess [|], |3j, that leads to spin flip scattering and, in semi- 
conductors without inversion symmetry, the Dyakonov- 
Perel (DP) process in which spin states precess be- 
cause of spin off-diagonal Hamiltonian matrix elements 
resulting from a combination of the spin-orbit interaction 
and inversion asymmetry. Typically, the DP mechanism 
dominates the spin dynamics in n-type III-V semiconduc- 
tors. An external magnetic field, in many cases required 
for control and manipulation in spintronics devices, can 
significantly influence electron spin dynamics. A mag- 
netic field has two main effects on electron spin relax- 
ation: (i) it quenches the DP process thereby tending 
to extend the spin lifetimes as a function of the mag- 
netic field 1^, and (ii) it opens an additional spin relax- 
ation process which tends to reduce the spin lifetimes 
in applied magnetic fields. |^ The latter process is due 
to the wave vector dependence of the conduction band 
(CB) electron g-factor. As a result of the variations in 
the (/-factor, electrons in different quantum states precess 
about a transverse magnetic field at different rates and 
thus lose spin coherence. For brevity we will refer to this 
process as a variable (7-factor (VG) mechanism. 

In contrast to previous studies ^, ^, ^, 0, ^, ||, 
p^ , pdl of spin relaxation in (bulk) n-type III-V semicon- 
ductors, we simultaneously treat the EY, DP, and VG 
processes on an equal footing and focus on the interplay 
between the various spin relaxation processes as a func- 
tion of the magnetic field. Thereby, we are able to study 
in detail the competition between the quenching of the 
DP process and the appearance of the VG process. 

Specifically, we calculate the longitudinal (Ti) and 
transverse (T2) spin relaxation times as a function of 



temperature, electron density and magnetic field. We 
find that the VG process dominantly infiuences the trans- 
verse (T2) spin relaxation time. In particular, as a result 
of the competition between the quenching of the DP pro- 
cess and the introduction of the VG process, there is a 
magnetic field for which the transverse (T2) spin lifetime 
is maximal. From the slope of T2 at small magnetic fields 
it is moreover possible to determine whether the DP or 
the EY process dominates spin relaxation at zero mag- 
netic field. In contrast, the magnetic field dependence of 
the longitudinal (Ti) lifetime is essentially unaffected by 
the VG process and dominated by the quenching of the 
DP process. Thus, they generally increase with field and 
saturate at a value given by the EY process. 

In an applied magnetic field, the CB electrons in a 
III-V semiconductor, e.g. GaAs, are described by the 
Hamiltonian 

where K — k — {e/hc)A{r) [A{f) is the vector po- 
tential], e{k) is the Kramers degenerate dispersion of 
CB electrons, hfl^ = fiBg*B is the Larmor frequency, 
h(ijA{K) = 25qk{K) is the splitting of the CB disper- 
sion due to the combination of spin-orbit interaction and 
inversion asymmetry, and 

nVLgiK) = 2aiK'^B + 2a^{K,B-K} + 2aQT{K, B) (2) 

is a term which gives rise to a wave vector dependence of 
the CB electron ^-factor. The definitions of the vectors 
K,{K) and t{K, B) and of the parameters 5q, a^, z = 4, 5, 6 
are given in Refs. |l^ and {, } indicates an anticom- 
mutator. 

Our calculation starts from the full quantum kinetic 
equations for the contour-ordered Green functions [Q, 
from which we derive, considering a classical homoge- 
neous magnetic field and using the fact that wave vector 
scattering is essentially instantaneous on the time scale 
of spin relaxation, a semiclassical kinetic equation for the 
CB electron density matrix. We then linearize this ki- 
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netic equation with respect to the CB electron spin den- 
sity, assuming, as an initial condition, small spin polar- 
ization. Treating scattering processes in the Born ap- 
proximation and expanding the collision integrals up to 
second order in the wave vector transfer (diffusion ap- 



proximation |1J]), we finally obtain a generalized Fokker- 
Planck-Landau equation for the spin density S{kt) which, 
in atomic units (with magnetic field along the z-axis) 
reads 



[dt - incCz ]S{kt) = pL + ^iA{k) + iigik)] 



with ilc the cyclotron frequency, Cz and the z- 
component and the squared total angular momentum op- 
erator in wave vector space, respectively, l/ri(A;) the sum 
of the (on-shell) wave vector relaxation rates for the vari- 
ous scattering processes, and V a differential operator in 
k — \k\ relevant to inelastic scattering processes. 

Equation (^) contains EY, DP, and VG processes and 
accounts for Larmor precession and orbital motion of the 
CB electrons in the magnetic field. More specifically, the 
EY process, due to genuine spin flip scattering events, is 
given by the tensor F, whereas the DP and VG processes 
originate from the interplay of spin conserving wave vec- 
tor scattering events described by the differential opera- 
tor T> — (1/2ti)£^ and the torque forces due to fljA and 
fig, respectively. The orbital motion encoded in —iflc^z 
and, to a lesser extend, the torque force due to i^L lead 
to a quenching of the DP process. 

It is possible to derive from Eq. (|^) general expressions 
for the spin relaxation rates without specifying whether 
the scattering processes are elastic {V = 0) or inelastic 
(2? 7^ 0). To that end, we follow Ref. [|j and employ a 
perturbative approach with respect to the torque forces. 
Our results are therefore valid for \iliA + ^g\Ti < 1. Ex- 
panding fljA{k) and flg{k) in terms of spherical harmon- 
ics Yi^ie, $), we find {i - 1, 2) H 

mr' ^ [Try' + [rrv' + [rrv" ^ (4) 

with the EY contributions (due to spin-fiip scattering) 
the DP contributions (due to inversion asymmetry) 

[TrV' = i\cI,\'fL + i\cl\'fL (6) 

[T,n-^= 2|C3^,pf|, + 2|C7|;pf|, + 2|C|'3pf|3, (7) 

and the VG contributions (due to the wave vector depen- 
dence of the CB electron g-factor) 

[Try' = 4\Dl\'r^„ (8) 



S(kt) + [!)- ——C^]S(kt) - r(k)S(kt), (3) 

2Tl(fc) 

I 

Here, Cl„^ and (i = X, Y, Z) are the expansion coeffi- 
cients of VL\^{k) and ^\{k), respectively, and {v = 1, 2, 3) 

= -Re(G(fc)Tr„(fc)), (10) 

where the brackets denote an average over k defined as 
((..._))= /„°°dfcfcV(fc)/»(-)/47r/„°°&fcV(fc)/(fc), 
with /(fc) = 1 — f{k) and /(fc) the equilibrium Fermi 
distribution function. The generalized wave vector re- 
laxation time Ti„-^(k) satisfies a differential equation 

[V + timilc + ^1) - ^Wrnik) = Clik), (11) 

with Ci{k) = CiAkHis + Cg{B)e[5i2 + Sio], ni - 

^l[Su1 - S,2], and 1/T;(fc) = l{l + l)/2Ti(fc). The con- 
stants characterizing the three spin relaxation processes 
are, respectively, Csf — 5'^{A + 2eg)i?oWo/2Aegm*, 
CiA = 26o/RoaQ, and Cg{B) — 2^,bB / R^aQ, where 
(5^ = 2A2/(A + eg){2A + Scg), A is the spin-orbit split- 
ting, €g is the band gap, Rq and Uq are the Rydberg 
energy and the Bohr radius, respectively, m* and mo are 
the CB electron mass and the mass of a bare electron, 
respectively, and /is is the Bohr magenton. The detailed 
form of V depends on the scattering processes and does 
not concern us here. ||l^ 

Note, as a consequence of the orthogonality of the angle 
dependences, the EY, DP, and VG spin relaxation rates 
are additive. The generalized relaxation rate l/T'/^(fc), 
on the other hand, is in general not proportional to the 
sum of the (on-shell) relaxation rates 1/T;(fc) because of 
inelasticity. A Matthiessen-type rule for 1/T;'^(fc) only 
holds for elastic scattering (see below). 

We are interested in the magnetic field dependence of 
the spin relaxation processes which, at least qualitatively, 
should not depend on the approximation adopted to de- 
cribe the scattering events. In the following, we treat 
therefore all scattering processes in the elastic approx- 
imation and neglect V in Eq. (pi]). Specifically, we 
take scattering on ionized impurities, acoustic phonons, 
and longitudinal optical (LO) phonons into account. The 
elastic approximation restricts our results to low enough 
temperatures, where electron-impurity scattering domi- 
nates, and to high enough temperatures, where electron- 
phonon scattering becomes essentially elastic. 
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FIG. 1: The top and bottom panels show, respectively, Ti 
and T2 in GaAs as a function of magnetic field for T = 
and n — lO^^cm"''. The contributions from the EY (long 
dash), DP (short dash) and VG (dot-dash) processes and the 
total relaxation time (solid) are shown in the main panel. The 
insets (same axis as the main panel) show the total relaxation 
times forn = 5 X lO^'^cm"^, 1 x lO^^cm"^, 5 x 10^** cm"^, and 
1 X lO^^cm"^ (top to bottom}. The squares and triangles are 
experimental data from Ref.[h9| at the respective densities. 



Within the elastic approximation Eq. (|ll| ) reduces to 
an algebraic equation which is readily solved to yield 

_ I ri{k) [C|^fc6j,3 + Cl{B)k\5n + Si,)] \ 

The k-average can be obtained either numerically or, at 
low and high temperatures, with saddle point techniques 
exploiting the peaked structure of the integrands. Within 
the elastic approximation it is sufficient to adopt the lat- 
ter. Details of the calculation will be given elsewhere. 

In Fig. |l|, we show calculated longitudinal (Ti) and 
transverse (T2) spin relaxation times for GaAs as a func- 
tion of magnetic field at T = and an electron den- 
sity of n = lO^^cm"'^. We show separately the contri- 
butions to the spin relaxation times from the EY, DP 



FIG. 2: The top and bottom panels show, respectively, Ti 
and T2 in GaAs as a function of magnetic field for T = lOOK 
and n — lO^^cm"^. The contributions from the EY (long 
dash), DP (short dash) and VG (dot-dash) processes and the 
total relaxation time (solid) are shown in the main panel. The 
insets (same axis as the main panel) show the total relaxation 
times for T = 150K, 200K, 250K, and 300K (top to bottom). 

and VG processes and the total spin relaxation time in- 
cluding all three spin relaxation processes. In the in- 
sets of Fig. ^ we give the total spin relaxation time 
for various electron densities at T = 0. The parameters 
needed to specify ^g{k) have been previously obtained 
partly experimentally by measuring combined cyclotron 
resonances (04, og) and partly theoretically within a five- 
level Kane model (a^): (04,05,06) = (97,-8,49) x 
10"24eVcm2Oe~\ ||T3| The parameter defining niA{k) 
is given by 60 = 0MhP/^{2m*feg. @ The remain- 
ing parameters, such as the effective CB electron mass 
or the deformation potential are available from standard 
data bases. Q 

For the temperature and density conditions in Fig. 
the electrons are degenerate and electron-ionized impu- 
rity scattering dominates. The VG process makes a small 
contribution to Ti which is dominated by the DP process 
at zero magnetic field. As the magnetic field is increased, 
the DP process is quenched. Thus, Ti increases monoton- 
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ically with increasing magnetic field saturating at high 
field at a value determined by the EY process which is 
not affected by the magnetic field. If the material pa- 
rameters had been such that the EY process dominated 
the DP process for Ti relaxation at zero magnetic field, 
Ti relaxation would not be significantly affected by the 
applied field. By contrast the VG process makes a sig- 
nificant contribution to T2 relaxation. At small applied 
magnetic fields the T2 lifetime increases with increasing 
magnetic field, but as the field continues to increase the 
VG process begins to dominate the relaxation so that T2 
has a maximum and begins to decrease for larger mag- 
netic fields. If the material parameters had been such 
that the EY process dominated the DP process for T2 re- 
laxation at zero magnetic field, the T2 relaxation would 
monotonically decrease with increasing magnetic field. 

The solid squares and triangles in the lower panel of 
Fig. are measured T2 spin lifetimes in GaAs at 5K 
from Ref. Ml at electron densities of 1 x lO^^cm"^ and 
5 X lO^^cni^. (Data for an electron concentration of 
1 X lO^^cm"-^ at 5K was also presented in Ref. |l^, but 
at this low density the electrons are bound to isolated 
donors and our theory does not apply.) There is good 
(order of magnitude) agreement between our calculation 
and these measured results, although there were no ad- 
justable parameters. Unfortunately, the magnetic fields 
in Ref. are not high enough to capture any effects 
due to the VG process. In particular, our prediction of 
the maximum of T2 remains to be experimentally veri- 
fied. 

In Fig. 1^, we show the various contributions to the 
Ti and T2 spin relaxation for GaAs as a function of 
magnetic field at T = lOOK and an electron density of 
n = lO^^cm"'^. In the insets of Fig. ||, we show the total 
spin relaxation time as a function of magnetic field for 
various temperatures at n = lO^^cm""^. For the temper- 
ature and density conditions in Fig. ^ the electrons are 
non-degenerate and electron-LO-phonon scattering is the 
dominant scattering process. As for the degenerate elec- 
tron case, the VG process makes a small contribution to 
Ti which is again dominated by the DP process at zero 



magnetic field. The DP process is quenched by the field 
so that Ti increases with field at small fields and satu- 
rates at a value determined by the EY process at large 
fields, jl^ Similar to the degenerate electron case, the 
VG process makes a substantial contribution to T2 re- 
laxation. At small fields the T2 lifetime increases with 
increasing field and at large fields the VG process begins 
to dominate the relaxation so that T2 has again a maxi- 
mum at some finite magnetic field. The sign of the slope 
in T2 at small magnetic fields is again a clear signature 
of whether the EY process (T2 decreases with increasing 
field) or DP process {T2 increases with increasing field) 
dominates T2 relaxation at zero magnetic field. Note, the 
qualitative behavior of the longitudinal and transverse 
spin relaxation times with increasing magnetic field is 
similar for degenerate and non-degenerate electrons, but 
the magnitude of the change is larger for non-degenerate 
electrons. 

In summary, based on a systematic kinetic approach, 
which treats the EY, DP, and VG processes on an equal 
footing, we calculated the longitudinal (Ti) and trans- 
verse {T2) spin relaxation times of CB electrons in n-type 
III-V semiconductors as a function of temperature, elec- 
tron density, and magnetic field. At finite magnetic field, 
the VG process competes with the DP and EY processes. 
We find that, as a consequence of the interplay of the DP 
and the VG processes, T2 can have a maximum as a func- 
tion of magnetic field. In contrast, Ti is not affected by 
the VG process and increases with magnetic field until it 
saturates at a value determined by the EY process. The 
sign of the change in T2 with increasing magnetic field at 
small fields indicates, moreover, whether the EY process 
or the DP process dominates T2 relaxation at zero mag- 
netic field. Our calculated results are in good agreement 
with existing experimental data in n-type GaAs and we 
make additional specific predictions for the magnetic field 
dependence of electron spin lifetimes that are subject to 
experimental check. 

This work was supported by the SPINs program of the 
U.S. Defense Advanced Research Projects Agency. 
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